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Abstract 

The HIJING and VENUS models of relativistic hadron-nucleus and nucleus- 
nucleus collisions are used to study interactions of proton and nuclei with nitro- 
gen, specific for the extensive air shower developments initiated by cosmic rays 
in the atmosphere. The transverse energytransverse momenta and secondary 
particles produced spectra as well as their energy and mass dependence were 
investigated in detail. Results are presented with particular emphasis on the 
contributions of minijets in HIJING model and validity of superposition models 
in this energy range . 

PACS : 13.85.Hd ; 13.85.Ni ; 13.85.Tp ;96.40.De;96.40.Pq;96.40.Tv 
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1 Introduction 



The investigation of the detailed shape of the energy spectrum and the mass composition of 
primary cosmic rays is currently a most active field of astrophysical research |lj . 

Experiments on satellites or with balloon-borne detectors give information up to energies 
of ca. 10 14 eV Due to their limitations in size and weight they can hardly be extended 
beyond 10 15 eV, and indirect techniques, the observation of the particle-cascades in the atmo- 
sphere (extensive air showers: EAS), have to be invoked. The information about nature and 
energy of the primary particles is reflected by the shower development || whose details and 
signatures for the primary particle depend on the high-energy nuclear interactions, governing 
the cascading processes. Thus the analysis requires a reliable description of these processes, 
formulated as a hadronic interaction model which can be used as generator of Monte-Carlo 
simulations of air showers. It should describe the currently available experimental informa- 
tion from accelerator experiments (in particular the data from the large collider facilities 
at CERN and Fermilab) and allow a justified extrapolation to experimentally unexplored 
energy regions. In the case of the EAS cascades, the quest is for the cross sections (multi- 
particle production, rapidity and transverse momentum distributuions) for hadron-hadron, 
hadron-nucleus and nucleus-nucleus collisions as function of the energy (from pion produc- 
tion threshold up to ultrahigh energies), most importantly for the forward fragmentation 
region, while actually the central region of the collisions is best studied at accelerators. The 
fragmentation region is nevertheless also relevant for the interaction models describing the 
experimental observations at SPS energies and beyond 

There are many hadronic transport models en vogue which address this problem. They 
comprise versions of the Dual-Parton model (DPM) [[F], Quark-Gluon String models (QGSM) 
g], and models designated with the name of the code like VENUS [|, FRITIOF 0, HIJING 



ll|]-[|T4| ,Parton Cascade Models [|T^],|T6[ and others. Some have been specifically developed 



as Monte-Carlo generator for air shower simulations at cosmic ray energies like DPM |17[ 
HEMAS and SYBILL fTJ. 

Recently [^(J the VENUS approach, linked to the CORSIKA code |^T] (now widely used 
for cosmic rays EAS simulations) has been used to scrutinize the superposition hypothesis 
in nucleus-nucleus collisions at ultrahigh energies. The hypothesis has been shown to be a 
rather good approximation, for the air shower cascade while Ranft |J using the DPMJET-II 
version of the DPM approach and considering the fragmentation and central regions with 
equal importance concluded that the superposition is a rather rough approximation of the 
reality. Our present work is based on the experience with a model and its extensions which 
are the basis of the HIJING code and used for an extrapolation of the particle production 
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dynamics from proton-proton (pp) to proton-nucleus (pA) and nucleus-nucleus (AA) inter- 
actions, taking into account the essential constraints of geometry and kinematics. At higher 
collision energies semi-hard processes are included by pertubative QCD processes (pQCD). 

While the HIJING model ignores final state interactions, the VENUS model includes rein- 
teractions of string segments among themselves and with spectator matter, with occurrence 
of " double strings" which mediate multinucleon interactions in nucleus-nucleus collisions. 

Both string-models have been applied to a variety of pp, pA and AA collision data (see 
references P|, ||12|| , |1l3| , [22]. However, a consistant intercomparison of predictions of multi- 



particle production, transverse momentum and rapidity distributions at ultrahigh energies 
and with respect to their relevance in the EAS cascade is missing. The present paper is a 
first attempt of such a comparison, revealing the most salient features and differences in a 
selected number of cases. We introduce the presentation of numerical results with a brief 
reminder of the basis of the HIJING model under consideration, stressing the different pro- 
cedures in defining the interacting nucleon configurations, the quark-gluon string formation 
and the decay into secondary particles. 

The models have been tested at accelerator energies for proton-proton and nucleus- 
nucleus interactions and then theoretical predictions on pseudorapidity distributions of trans- 
verse energytransverse momenta and secondary particles spectra as well as their energy and 
mass dependence are given using HIJING model for proton - Air Nucleus (p+Air) interactions 
between 1 TeV - 1000 TeV and for Nucleus - Air (A+Air) interactions at 17.86 TeV/Nucleon 
corresponding to 1 PeV for iron (Fe) nucleus. A comparison with recent results |2(| at the 
same energy using VENUS model is also done. We have investigated Feynman scaling be- 
haviour of the model in this energy region ,and the multiple minijets production is considered 
for study of the charged multiplicity distributions. Finally a brief discussion on validity of 
superposition models ,taken into consideration mean integrated values of transverse energy 
predicted by HIJING model, is presented. 

2 Outline of HIJING Model 



A detailed discussion of the HIJING Monte Carlo model was reported in references |l 1 
0. The formulation of HIJING was guided by the LUND-FRITIOF and Dual Parton 
Model(DPM) phenomenology for soft nucleus-nucleus reactions at intermediate energies 
(a/s < 20 GeV) and implementation pQCD processes in the PHYTHIA model[|23] for 
hadronic interactions. We give in this section a brief review of the aspect of the model 
relevant to hadronic interaction: 
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1. Exact diffuse nuclear geometry is used to calculate the impact parameter dependence 
of the number of inelastic processes. 



Soft beam jets are modeled by quark-diquark strings with gluon kinks along the lines 
of the DPM and FRITIOF models. Multiple low px exchanges among the end point 
constituents are included. 

The model includes multiple mini-jet production with initial and final state radiation 
along the lines of the PYTHIA model and with cross sections calculated within the 
eikonal formalism. 



4. The hadronization of single chains is handled by the LUND code JETSET 7.3 [f23|1 ,that 
summarizes data on e + e~. In this picture , each nucleon-nucleon collision results in 
excitation of the nucleon by the stretching of a string between the valence quark and 
diquark ( longitudinal excitation , without colour exchange). A phenomenological exci- 
tation function determines the mass and momentum of the string after each interac- 
tion. After the last interaction the string decays to produce particles. 

5. HIJING does not incorporate any mechanism for final state interactions among low 
Pt produced particles nor does it have colour rope formation. 

Unlike heavy ion collisions at the existing AGS/BNL and SPS/CERN energies, most of 
the physical processes occurring at very early times in the violent collisions of heavy nuclei at 
cosmic ray energies involve hard or semihard parton scatterings which will result in enormous 
amount of jet production and can be described in terms of pQCD. Assuming independent 
production, it has been shown that the multiple minijets production is important in proton- 
antiproton (pp) interactions to account for the increase of total cross section and the 
violation of Koba-Nielsen-Olesen (KNO) scaling of the charged multiplicity distributions 

0- 

In high energy heavy ion collisions, minijets have been estimated [^J to produce 50% 
(80%) of the transverse energy in central heavy ion collisions at RHIC (LHC) energies. 
While not resolvable as distinct jets, they would lead to a wide variety of correlations, as 
in proton-proton (pp) or antiproton-proton pp collisions, among observables such as multi- 
plicity, transverse momentum, strangeness, and fluctuations that compete with the expected 
signatures of a QGP. Therefore, it is especially important to calculate these background 
processes. 
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These processes are calculated in pQCD starting with the cross section of hard parton 
scatterings pBj : 



dfT ' ! ' 1 - AT £ xix 2 f a ( Xl , P|)/ 6 (x 2 , P 2 )da ab (s, t, u)/dt, 



a.b 



where the summation runs over all parton species, 1/1,1/2 are the rapidities of the scattered 
partons and xi,x 2 are the fractions of momentum carried by the initial partons and they 
are related by x x = x T {e m + e m )/2, x 2 = x T {e~ m + e~ m ), x T = 2P T /y/s. A factor, K « 2 
accounts roughly for the higher order corrections. The default structure functions, f a {x, Q 2 ), 
in HIJING are taken to be Duke-Owens structure function set 1 P7| . 

Integrating Eq. |I| with a low Pt cutoff Pq, one abtain the total inclusive jet cross section 
a j e t. The average number of semihard parton collisions for a nucleon-nucleon collision at 
impact parameter b is aj et T^(b), where T^(b) is partonic overlap function between the two 
nucleons. In terms of a semiclassical probabilistic model |?4], |12| , the probability for multiple 
minijets production is : 

g 3 (b) = h^Sl e -^ N{b)j ■ > L (2) 

Similarly, the soft interactions are represented by an inclusive cross section a so f t which, unlike 
a j e t, can only be determined phenomenologically. The probability for only soft interactions 
without any hard processes is : 

g (b) = [1 - e" T ^ TN( f ,) ]e-^ Mb) . (3) 

and the total inelastic cross section for nucleon-nucleon collisions : 



/oo 

j=0 

= J d 2 b[l - e -^°t t+r7 ^ )TN{b) }. (4) 
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Define a real eikonal function, 



X(b, s) = ^a soft (s)T N (b, s) + ^a jet (s)T N (b, s), 



(5) 



the elastic, inelastic, and total cross sections of nucleon-nucleon collisions are given by : 



a,i = it db 



1 -e 



-x(M) 



(6) 



(Ti n = IT db 2 
JO 

POD 

a = 2n db 2 
Jo 



1 _ e -2x(M) 
1 - e~ x(M; 



(7) 



Assuming that the parton density in a nucleon can be approximated by the Fourier transform 
of a dipole form factor, the overlap function is : 



T N (b, s) 



, Xo(fl 

'(?soft{s)' 



(9) 



with 



96 



(10) 



where /io = 3.9 and vt6q(s) = (t = cr S oft(s)/2 is a measure of the geometrical size of the 
nucleon. The eikonal function can be written as : 



X(b, s) = x(£, s) = [1 + a jet (s)/a soft (s)]xo(0- 



(11) 



Po — 2 GeV/c and a constant value of cr so ft(s) = 57 mb are chosen to fit the experi- 
mental data on cross sections [Q in pp and pp collisions. The equations listed above are 
used to simulate multiple jets production at the level of nucleon-nucleon collisions in HI- 
JING Monte Carlo program. Once the number of hard scatterings is determined, PYTHIA 
algorithms generate the kinetic variables of the scattered partons and the initial and final 
state radiations. 

After all binary collisions are processed, the scattered partons in the associated nucleons 
are connected with the corresponding valence quarks and diquarks to form string systems. 
The strings are then fragmented into particles. 



6 



The HIJING model incorporate nuclear effects such as parton shadowing and jet quench- 
ing. HIJING is designed also to explore the range of possible initial conditions that may occur 
in nuclear collisions at Cosmic Ray and colliders (RHIC,LHC) energies. 

To include the nuclear effects on jet production and fragmentation, the EMC [28|] effect 



of the parton structure functions in nuclei and the interaction of the produced jets with the 
excited nuclear matter in heavy ion collisions are considered. 

It has been observed ]28| that the effective number of quarks and antiquarks in a nucleus 
is depleted in the low region of x. While theoretically there may be differences between 
quark and gluon shadowing,one assume that the shadowing effects for gluons and quarks are 
the same. 

At this stage, the experimental data unfortunately can not fully determine the A de- 
pendence of the shadowing. In the HIJING model the A dependence are taken from Ref. 



Ra(x) = 



fa/A(x) 



Af a/N {x) 

= 1 + 1.19 ln 1/6 A [a; 3 - 1.5(x + x L )x 2 + 3x x L x] 
. 1.08(^/3-1) 

- [aA ~ HATY) v " (12) 

a A = 0.104 1 / 3 -!), (13) 



where Xq = 0.1 and xl = 0.7. The term proportional to a a in Eq. [12] determines the 
shadowing for x < Xq with the most important nuclear dependence, while the rest gives the 
overall nuclear effect on the structure function in x > Xo with some very slow A dependence. 

To take into account impact parameter dependence, one assume that the shadowing effect 
aA is proportional to the longitudinal dimension of the nucleus along the straight trajectory 
of the interacting nucleons. The values of in Eq. [12] are parametrized like : 

4 



a A (r) = 0.1(^/3 _ l)-^/l- r yRl (14) 

where r is the transverse distance of the interacting nucleon from its nucleus center and 
Ra is the radius of the nucleus. For a sharp sphere nucleus with overlap function T^(r) = 
{3A/2ttR a )Ji -r 2 /R% the averaged a A (r) is a A = n ff A dr 2 T A {r)a A {r)/A. 



7 



To simplify the calculation during the Monte Carlo simulation, one can decompose 
Ra(x, r) into two parts, 

R A (x, r) = R° A (x) - a A (r)R s A (x), (15) 

where a A (r)R s A (x) is the term proportional to a / i(r) in Eq. [12] with a A (r) given in Eq. p~4] 
and R A (x) is the rest of Ra(%, r). Both R A (x) and -R^(x) are now independent of r. 

The effective jet production cross section of a binary nucleon-nucleon interaction in A+B 
nuclear collisions is then, 



vfJ(rA, r B ) = Vjet - a A (r A )af et - a B (r B )af et + a A (r A )a B (r B )af e f } (16) 



where cr? et , a^ et , a^ et and cr^f can be calculated through Eq. |T] by multiplying 
f a (xi, P|)/f,(x 2 , P|) in the integrand with R A (xi)R° B (x 2 ), R A (x 1 )R° B (x 2 ), R A (x 1 )R s B (x 2 ) and 
R s A (xi)R s B (x2) respectively. With calculated values of a° et , af et) af et and afjj* , one get the 
effective jet cross section cr|f/ for any binary nucleon-nucleon collision. 

Another important nuclear effect on the jet production in heavy ion collisions is the final 
state integration. In high energy heavy ion collisions, a dense hadronic or partonic matter 
must be produced in the central region. Because this matter can extend over a transverse 
dimension of at least R A , jets with large Pt from hard scatterings have to traverse this hot 
environment. For the purpose of studying the property of the dense matter created during 
the nucleus-nucleus collisions, HIJING model include an option to model jet quenching in 



terms of a simple gluon splitting mechanism |T2"f,||L3||. 



The induced radiation in HIJING model is given via a simple collinear gluon splitting 
scheme with energy loss dE/dz. The energy loss for gluon jets is twice that of quark jets||30||. 
One assume that interaction only occur with the locally comoving matter in the transverse 
direction. The interaction points are determined via a probability : 

rfP= ^ e -^/A % (17) 

with given mean free path A s , where I is the distance the jet has traveled after its last 
interaction. The induced radiation is simulated by transferring a part of the jet energy 
AE(£) = £dE/z as a gluon kink to the other string which the jet interacts with. The 
procedure is repeated until the jet is out of the whole excited system or when the jet energy 
is smaller than a cutoff below which a jet can not loss energy any more. This cutoff is 
taken as the same as the cutoff Pq for jet production. To determine how many and which 
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excited strings could interact with the jet, one assume a cross section of jet interaction so 
that excited strings within a cylinder of radius r s along the jet direction could interact with 
the jet. X s should be related to r s via the density of the system of excited strings. The values 
for X s and r s are taken as two parameters in the model. 

The main usefulness of these schematic approaches for nuclear shadowing and jet quench- 
ing is to test the sensitivity of the final particle spectra . 

For a detailed description of the VENUS model as applied in the calculations presented 
in this paper we refer to excelent review of Werner (sections 6-10) §J. 



3 NUMERICAL RESULTS 

For used a Monte Carlo generator in Cosmic Ray physics this model should provide the basic 
hadronic interaction term for the Cosmic Ray cascade i.e should provide the cross sections for 
hadron-hadron,hadron-nucleus and nucleus-nucleus collisions as function of the energy . Also 
should provide a good description of secondary particles production since, secondary 7r° and 
rj mesons are the source of the electromagnetic shower , secondary and K mesons are the 
source of Cosmic Ray Muons and the source of atmospheric Neutrinos produced by Cosmic 
Ray cascade and secondary charmed mesons are the source for prompt Muons and Neutrinos. 
The model should work from accelerator energies up to the highest possible primary energies . 
Experimental observations like rapidity plateaus and average transverse momenta rising with 
energyKNO scaling violation, transverse momentum - multiplicity correlations and minijets 
pointed out that soft hard and processes are closely related and all these properties were 



understood within the HIJING model 11-14 



The HIJING model provides a framework not only for the study of hadron-hadron inter- 
actions ,but also for the description of particle production in hadron nucleus and nucleus - 
nucleus collisions at ultrahigh energies. 

The relevance of an event generator like HIJING or VENUS for hadron production cross 
sections in the Cosmic Ray energy region can only be claimed if the model agrees to the 
best available data in the accelerator energy range and if it shows a smooth behaviour in the 
extrapolation to ultrahigher energies . 

For the Cosmic Ray Cascade in the atmosphere only hadron - nucleus and nucleus- 
nucleus collisions are relevant with Nitrogen 14 N beeing the most important target nucleus. 
However experimental data are of much better quality in hadron-hadron, and especially in 
proton-proton collisions or antiproton-proton collisions than for collisions of hadrons with 
light nuclei. The scarcity of the data in this region was also pointed out in studing strangeness 
production at SPS energies 
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So we start with the study of proton - proton collisions and nucleus - nucleus collisions 
at accelerator energies . 



3.1 PROTON - PROTON AND NUCLEUS-NUCLEUS INTER- 
ACTIONS AT ACCELERATOR ENERGIES 

We used the program HIJING with default parameters: 

1. IHPR2(11) = 1 gives the baryon production model with diquark-antidiquark pair pro- 
duction allowed, initial diquark treated as unit; 

2. IHPR2(12)=1, decay of particle such as tt° , K® , A , £ , S , £1 are allowed ; 

3. IHPR2(17)=1 - Gaussian distribution of transverse momentum of the sea quarks ; 

4. IHPR2(8)=0 - jet production turned off for theoretical predictions denoted by HIJING 
model ; 

5. IHPR2(8)=10-the maximum number of jet production per nucleon-nucleon interaction 
for for theoretical predictions denoted by HIJING^ at SPS energies and 300 GeV ; 

6. IHPR2(8)=20 for energies > 1 TeV . We have neglected also nuclear shadowing ef- 
fect, and jet quenching (see section 2) and we have a cut in pseudorapidity r\ > 0.5 for 
ultrahigh energies (> 1 TeV). 

In Table I the calculated average multiplicities of particle at Ei^ = 200 GeV in (pp) in- 
teractions are compared to data. The theoretical values as predicted by the model HIJING 
and HIJING^ are obtained for 10 5 generated events and in a full phase space. The val- 
ues HIJING^ include the very small possibility of mini jet production at these low SPS 
energies. The experimental data are taken from Gazdzicki and Hansen . The theoretical 
values given by VENUS model(as computed here) are obtained for 10 4 generated events 
and in a full phase space. The values labeled DPMJET II are taken from Ranft ||. 

The small kaon to pion ratio is due to the suppressed strangeness production basic to 
string fragmentation. Positive pions and kaons are more abundant than the negative ones 
due to charge conservation. We note that the integrated multiplicities for neutral strange 
particle < A >, < A >, < > are reproduced at the level of three standard deviations for 
pp interactions at 2QQGeV and 300 GeV. (see also Fig. la and Fig. lb -experimental data 



from Lo Pinto et al. [33] for pp interactions at 300 GeV.) However the values for < p > 



and < A > are significantly over predicted by the models. 
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For completeness we include a comparison of hadron yields at collider energies Elab = 
160 TeV for pp interactions, where mini-jet production plays a much more important role. 
From different collider experiments Alner et al. (UA5 Collaboration) []33| attempted to 



piece together a picture of the composition of a typical soft event at the Fermilab SppS 
collider ||35|| . The measurements were made in various different kinematic regions and have 
been extrapolated in the full transverse momentary) and rapidity range for comparison 
as described in reference ||34|| . The experimental data are compared to theoretical values 



obtained with HIJING^ in Table II. It was stressed by Ward |35[ that the data show a 
substantial excess of photons compared to the mean value for pions < tt + + tt~ > . It was 
suggested as a possible explanation of such enhancement a gluon Cerenkov radiation emission 
in hadronic collision [|36[] . Our calculations rules out such hypothesis. Taking into account 
decay from resonances and direct gamma production, good agreement is found within the 
experimental errors. 

In the following plots the kinematic variable used to describe single particle properties 
are the transverse momentum px and the rapidity y defined as usual as: 

1 , E + p 3 . E + p 3 . . 

y = —In = In (18) 

y 2 E-p 3 m T K J 

with E,p 3 , and being energy, longitudinal momentum and transverse mass = 



m o + Pt w hh mo being the particle rest mass. 
The pseudorapidity rj is used rather than the rapidity since for r\ no knowledge of particle 
masses is required. 

77 = -In - ^ 3 = —Intan- (19) 
' 2 p-p 3 2 1 ; 

where p is the projectile nucleon momentum and 9 is the scattering angle. 
Feynman xp variable is defined for ultrahigh energy as: 

x F = 2^fsinh{y cm ) (20) 
V s 

where y cm ans s are rapidity and total energy in center of mass frame (cms). 

In Figure 1 we compare rapidity and transverse momentum distributions for strange 
particles in proton-proton interactions at 300 GeV given by HIJING model with experimental 
data |33|]. The agreement is quite good. However ,it will be interesting to investigate in the 
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future the Feynman scaling behaviour of the model at the accelerator energies since the 
forward fragmentation region seems to play an important role ||, 0,0. 

This year marked a real progress in the field of high energy nuclear collisions .Now 
an entirely new domain of energies has become accessible with heavy nuclear beams [37 |. 
Recently J22| using HIJING and VENUS models the systematics of strangeness enhancement 
was calculated for pp,pA and AA interactions at SPS energies and it was stressed out that 
the enhancement of strangeness has its origins in non-equilibrium dynamics of few nucleons 
systems. To clarify the new physics much better quality data on elementary pp as well as 
other light ion reactions will be needed for tunned models . 

We foccused our analyses concerning particle production in nucleus - nucleus interactions 
at SPS energies , where the models should be better tested , mainly for two simmetrical 
interactions S + S , Pb + Pb and for asymmetrical one S + Pb. 

It was shown that the negative pion rapidity densities are well accounted for both HIJING 
and VENUS @j, but the flat valence proton distribution in S + S is only reproduced by 
VENUS. Recall that VENUS includes a model of final state interactions in dense matter as 
well as a colour rope effect , called double strings. VENUS predicts a much higher degree of 
baryon stopping at midrapidity than HIJING . These conclusions are confirmed also for our 
calculations in S + S at 200 AGeV (see Figs. 2a, 2b) and Pb + Pb collisions at 160 AGeV 
(see Figs. 2c, 2d) for A and A . The experimental data are taken from Alber et al. p8| . 
The rapidity distributions for antiproton p and for negative kaons K~ are represented for 
S + S interactions (Fig.3a(p) and Fig.3b(i^ - )) and Pb + Pb interactions (Fig.3c(p) and 
Fig. 3d {K~)) in comparison with some experimental data taken from Baechler et al. |39| for 



K~ and from Murray et al. [4(J for p . For asymmetrical interactions S + Pb we predict in 
both models the rapidity spectra for 7r~ (Fig. 4a) , K~ (Fig. 4b) , A (Fig. 4c) and A (Fig. 
4d )) at 200 AGeV . We mention that we have used VENUS model with the option without 
decaying resonances (dotted histograms). If will allow decay of resonances then we get an 
increase of the values especially at mid rapidity which improve agreement between theory 
and experimental data (dashed histograms). 

It will be interesting to compare these results with upcoming data to test if the strangeness 
enhancement increases from SS to SPb or from SS to PbPb interactions. 



3.2 PROTON -AIR NUCLEUS INTERACTIONS AT ULTRA- 
HIGH ENERGIES 

Since not enough data are avaiable in the fragmentation region of hadron collisions with 
light target nuclei,many features of particles production in collisions involving nuclei can 
only be extracted from the study of models. These kind of analysis have been done recently 
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in VENUS model |2(J and DPM model In HIJING model some specific interactions 



were investigated at RICH and LHC energies ||T2f-||14||. 



Taking into account the traditional concept of collisions between nucleons and air nuclei 
the primary loses its energy during the interaction and this fluctuate between zero and 
100 % . The fraction of the energy lost and used for producing new particles is then referred 
to as inelasticity. We do not investigate traditionaly inelasticities (baryonic, electromagnetic 
or mesonic) see |H| . 

We generate for proton-Air Nucleus (p + Air) interactions 10 4 minimum bias events 
{brnin = 0, b max = 5fm) . In order to give an idea about energy used for producing new 
particle we investigate in Figure 5 the transverse energy pseudorapidities spectra and their 
dependence with energy for all secondaries (Fig. 5a), all neutral (Figure 5b) ,all charged 
(Fig. 5c) and gluons (Fig. 5d). As we see from Figure 5 gluons carried an important fraction 
from transverse energy and this fraction increase with increasing energy (from 2.8 % at 1 
TeV to 17.3 % at 1000 TeV). Also the percentage of occurence of gluons increase from 7.20 
% at 17.86 TeV to 10.71 % at 1000 TeV (see table 3). 

The rapidity distributions and energy dependence of main secondary particles are shown 
in Figure 6 and Figure 7. Secondary (Fig. 6a and 6b) and (Fig. 6c and 6d) are the source 
of Cosmic Ray Muons and the source of atmospheric Neutrinos produced by the Cosmic 
Ray cascade. Secondary 7r°(Fig.7a) and rj mesons are the source of the electromagnetic 
shower and secondary charmed mesons (Fig. 7d) are the source for prompt Muons and 
Neutrinos . The lambdas have been shown (Fig. 7c), because they can be produced from 
protons by exchanging a single valence quark by strange quark. The outer maximum of their 
distributions are due to this process. This component is identified in all nucleon distributions 
(neutron,proton). The statistics (10 4 events generated) seems not to be enough for charmed 
mesons (D^) but we give such distributions only to show that introducing HIJING code 
in a shower code with a much higher statistics at simulation level a study of promt muon 
component should be feasible. 

We see from Figure 6 that mesons distributions exhibit a broad structureless shape which 
does not depend strongly on the type of mesons but have a dependence on energy from 
1 TeV to 1000 TeV which is more pronunciated compared with VENUS results (see reference 

More less energy dependence is seen in tranverse momenta distributions for secondaries 
in p+Air Nucleus interactions in Figure 8( for all charged particles - Fig. 8a, for proton - 
Fig. 8b, for positive pions 7r + - Fig. 8c and for negative pions 7r~ -Fig.8d) 

Trying to stress out the relevance of accelerator data on particle production in hadron 
- nucleus collisions for Cosmic Ray cascade we study the Feynman scaling behaviour of 
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p + Air — > tt^ + X and p + Air — > p + p + X in Fig. 9a and Fig. 9b respectively. We plot 
the xfcLN/cLxf distributions for laboratory energies of 1 TeV (dotted histograms), 100 TeV 
(dashed histograms), 1000 TeV (solid histograms). The violations of Feynman scaling which 
occur are connected with known rise of rapidity plateau for all kinds of produced particles 
and with production of minijets. Due to minijets Feynman scaling is more strongly violeted 
especialy in the region xf > . The violation of Feynman scaling are less dramatic in 
DPMJET II model ||,|| and appear only arround xf = and xf = l.We note also that 
HIJING show violation of KNO scaling due to the production of multiple minijets and the 
tendency becomes stronger with increasing energy |12| . 



In order to evaluate multiplicity distributions for the charged particles in p + Air inter- 
actions (Fig. 10a) we differentiate the contributions from soft (Fig. 10b events with Nj et = 0) 
and hard processes (Fig. 10c events with Nj et = 1 and Fig.lOd events with Nj et > 1 ) where 
Nj et is the number of minijets produced in that events. Our calculations including the ef- 
fects of multiple minijets are the contributions from the events which have hard colisions 
with Px > Pq- Analysing Figure 10 it is clear that the events at the tails of the charged 
multiplicity distributions in p+Air interactions are mainly those with multiple minijets pro- 
duction. 

3.3 NUCLEUS-AIR NUCLEUS INTERACTIONS AT ULTRA- 
HIGH ENERGIES 

Due to the large fraction of nuclei in primary Cosmic Rays , Nucleus-Air collisions are of 
great importance in the EAS development.lt is important ,that the model will be able to give 
a good description of hadron production in nucleus-nucleus interactions. So, in this subsection 
we try to investigate mainly the dependence on projectile mass for specific interactions for 
EAS (He,Ne,S,Fe+Air Nucleus) at 17.86 TeV/Nucleon which correspond to 1 PeV laboratory 
energy for Fe nucleus . 

For a real comparison ,the nucleus-nucleus collisions geometry should be the same for all 
of the Monte Carlo codes since nuclear density distributions are well known from nuclear 
physics Hi] .At SPS energies the calculated number of target and projectile participants as 



well as the number of participants as a function of reaction impact parameter (b fm) shows 
no difference in HIJING and VENUS model §. 

The results depicted in Figure 11 are obtained for 10 4 generated events and for the follow- 
ing intervals {b m i n —b max ) of impact parameter : Fe+Air (0-13 fm), S+Air(0-ll fm),Ne+Air(0- 
10 fm), He+Air (0-7 fm) - dotted histograms ; Fe+Air (0- 8 fm) ,S+Air(0-7 fm),Ne+Air(0- 
6 fm), He+Air (0-5 fm ) -dashed histograms ; Fe+Air (0-5 fm),S+Air(0 - 5 fm),Ne+Air(0 - 
5 fm), He+Air (0-4 fm) -solid histograms. 
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From Figure ll(a,b,c,d) we see that at 17.86 TeV/Nucleon the results have strongly 
dependence on impact parameter intervals (b m i n -b max ) and also the possibility of mixing 
interactions if we are interested in rapidity spectra only. 

Therefore for comparison of HIJING and VENUS results at 17.86 TeV/Nucleon for 
p+Air, He+Air, Ne+Air,S+Air,Fe+Air and for p+Air at 1000 TeV we try to get approxi- 
matively the same number of participants. For HIJING model we get the values for mean 
numbers of participants listed in last lines of Table 3 for 10 4 generated events and the follow- 
ing impact parameter intervals (b m i n — b max ): He+Air(0-4 fm),Ne+Air(0-6 fm),S+Air (0-7 
fm),Fe+Air (0-8 fm). p + Air (0- 5 fm). The results for VENUS model are taken from Schatz 



et al. 20 



Table 3 lists the frequency of various particles among the secondaries for collisions of 
different nuclei with nitrogen. Our calculations confirm the results from reference |20 .The 
percentages do not seem to depend strongly on projectile mass nor, as shown by proton results 
on energy. We see only a slight tendency of increasing strangeness production with increasing 
primary mass for VENUS results. Analysing the results of table 3 we see that the VENUS 
model predict more pions and kaons , but less gamma particles that HIJING model . We 
remark also differences in total multiplicities at ultrahigh energy which can not be explained 
only by the difference between total number of participants . It seems that "double string" 
mechanism change considerably the baryon spectra and allows the baryon number to migrate 
several units of rapidity from the end point rapidity. 

Important for Cosmic Ray studies and EAS development is also the dependence of par- 
ticle production on the nuclear target and projectile. Taken into consideration the same 
conditions as those reported for the values listed in Table 3 the theoretical predictions for 
mass dependence of pseudorapidities spectra are given for transverse energy in Figure 12, for 
transverse momenta in Figure 13 and for main secondaries produced new particles in Figure 
14. The transverse momenta distributions of all charged, proton, 7r + ,7r _ are represented in 
Figure 15. For Feynman xpdN/dxF distributions we give theoretical values for specific EAS 
interactions in Figure 16 only for charged pions and for sum of protons and antiprotons . We 
see a slight mass dependence of Feynman distributions and transverse momenta distributions 
for A>20. 

Instead of the proper sampling of Nucleus + Air Nucleus scattering events , an approxi- 
mation often applied in EAS development is the so called superposition model. There are two 
different possible superposition models: a nucleus-nucleus collision A-B with N part partici- 
pating nucleons is approximated as the superposition of N part simultaneous nucleon-nucleon 
collisions and the second one : a nucleus-nucleus collision A-B with N pro j participating 
projectile nucleons is approximated as the superposition of N pro j simultaneous nucleon-B 
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collisions. The validity of this principle was analysed in some recent works ||20 |, ||,|J with 
different conclusions. 

We investigate in HIJING model the integrated mean transverse energy for secondaries 
produced in Nucleus+Air Nucleus interactions at 17.86 TeV/Nucleon. So, we generate 10 4 
events in the same impact parameter interval (0-5 fm) for Fe,S,Ne,He,p+Air interactions. 
The values obtained for mean projectile participants (N pro j) and mean number of binary colli- 
sions N co u (which include nucleon-nucleon(iV— N), nucleon-wounded nucleon(iV— N W ),N W — N 
and N w — N w collisions) are listed in Table 4. The results given in Table 4 are for all sec- 
ondaries , all charged and all neutrals particles produced in interactions. The values of 
integrated mean transverse energy < Et > scale with number of binary collisions in nucleus- 
nucleus interactions at this energy N co u and scaling proprietes are valid for A > 20 and should 
not be applied to ligth nucleus+Air interactions. Since integrated mean value of transverse 
energy < Et > is a measurable quantity it will be interesting to verify this scaling at RHIC 
energies. 

In order to evaluate multiplicitity distributions for charged particles in A+Air interactions 
we differentiate the contribution for all events (soft + hard) (Fig.l7a for Fe+Air and Fig. 17b 
for S+Air ) and only hard processes - events with Nj et > 1 (Fig. 17c for Fe+Air and Fig.l7d 
for S+Air ). We can see from Figure 17 that the low multiplicity events are dominated by 
those of no jet production while high multiplicity events are dominated by those of at least 
one jet production. Also it is clear from Figure 17 that the contributions from the events 
which have hard collisions increase with increasing of available energy and of projectile mass. 

At ultrahigh energies nuclear effects like nuclear shadowing of partons and jet quench- 
ing (see section 2) , should have important contributions [0],[13|. In HIJING model 



a simple parametrisation of gluon shadowing (see eqs 12-15) and a schematic quenching 
model (see eqs 16-17) were introduced to test the sensitivity of the final distributions to 
these aspects of nuclear dynamics. Trying to estimate the results taken into consideration 
and neglected such effects we expect bigest differences for more central Fe+Air interac- 
tions at 17.86 TeV/Nucleon. Therefore we will estimate these effects for impact parameter 
interval(0 — 5/m)). In Figure 18 we give theoretical values for pseudorapidity distribu- 
tions for charged pions (Fig. 18a and Fig. 18b) and for charged kaons (Fig. 18c and Fig.l8d) 
without parton shadowing and jet quenching (solid histograms) and with parton shadow- 
ing and jet quenching(dashed histograms). The quenching mechanism in HIJING is limited 
to pt > 2 GeV/c partons. For pseudorapidity distributions of main secondary particles 
produced in EAS specific interactions in this energy region these effects should be neglected. 

However ,these effects are more pronunciated at the early stage of collisions Jike we see 
from Figure 19 which depicted pseudorapidity distributions for gluons (Fig. 19a) and quark- 
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antiquark pairs(Fig.l9b) as well as transverse momenta distributions for valence and sea 
partons (Fig. 19c) and valence parton only (Fig.l9d). Of great importance are also energy 
densities which are produced in these interactions ,at the early collisions time but their values 
are under study now. 



4 Conclusions 

In this paper we have performed an analysis of particle production in pp and AA colli- 
sions at SPS CERN-energies using the HIJING and VENUS models and investigate possible 
applications of the models at Cosmic Ray energies mainly for p+Air Nucleus at 1 TeV 
,10 TeV, 1000 TeV and Nucleus-Air Nucleus interactions (Fe,S,Ne,He+Air,p+Air) at 17.86 
TeV/Nucleon . 

We stressed out that integrated multiplicities are quite well reproduced at the level of 
three standard deviations for pp interaction at 200 GeV and 300 GeV. However the values 
for < p > and < A > are significantly overpredicted by the models. A very good agreement 
is found within experimental errors for ultrahigh energies ( 160 TeV) in HIJING approach 
where mini-jet production plays a much more important role. 

For nucleus-nucleus collisions at SPS energies the rapidity spectra are well accounted for 
both HIJING and VENUS models for mesons ,but VENUS model seems to give a more good 
description for strangeness production. At SPS energies a final state interactions in dense 
matter as well as a colour rope effect predicts much higher degree of baryon stopping at 
midrapidity than HIJING. 

The event generator VENUS |[2(| and HIJING in the present study were tested to simulate 



ultrahigh energy collisions specific for EAS development. The transverse energytransverse 
momenta and secondary particles produced spectra as well as their energy and mass depen- 
dences were investigated in detail. 

The contributions from the events which have hard collisions are stongly evidentiated at 
the tail of charged particles multiplicity distributions and increase with increasing energy 
and increasing mass of projectile. Feynman scaling is strongly violeted in HIJING model 
due to multiple minijets events. For Cosmic Ray energy region and specific EAS interactions 
the effects of parton shadowing and jet quenching should be neglected. 

Possible utilization inside a shower code of these models allows to extend the analyses 
and to release the simple superposition model which is not valid at least for integrated mean 
transverse energy of secondaries particle . The thoretical values predicted by HIJING model 
suggest a scaling with number of binary nucleon-nucleon collisions. It will be interesting to 
verify this hypothesis in the future experiments at RHIC energies. 
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Figure Captions 



Figure 1: Rapidity (Fig. la) and transverse momenta(Fig.lb) distributions of A (dotted 
histograms), A (dashed histograms), and Kg (full histograms) produced in pp interactions 
at 300 GeV . HIJING results are shown by histograms. The experimental data are taken 
from p3|. 



Figure 2: Rapidity distributions of A and A produced in central SS collisions at 200 AGeV 
(Fig. 2a and Fig. 2b) respectively and in central PbPb collisions at 160 AGeV (Fig.2c and 
Fig. 2d respectively). Expectations based on HIJING model are depicted as solid histograms 
. The theoretical predictions based on VENUS model are depicted as dotted histograms 
(option without decaying of resonances) and as dashed histograms (option with decay of 
resonances). The NA35 data (full circles) are from Alber et al. J5H]. The open circles show 
the distributions for SS collisions reflected at yi a b = 3.0 . 
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Figure 3: Comparison of central S + S collisions at 200 AGeV (Fig. 3a -antiproton and 
Fig. 3b - negative kaons) with central PbPb collisions at 160 AGeV (Fig.3c-antiproton and 
Fig. 3d- negative kaons). The experimental data are from NA44 [4~0 for antiprotons and from 



3£| for negative kaons. The open circles show the distributions for SS collisions reflected 
at yiab — 3.0 . The solid, dashed and dotted histograms have the same meaning as in Figure 
2. 



Figure 4: The various theoretical predictions for S + Pb — ► 7r + A (Fig. 4a), S + Pb — ► 
K~ + A(Fig.4b), S + Pb -> A + X (Fig.4c) and S + -> A + X (Fig.4d) at 200 AGeV. 
The solid and dotted histograms have the same meaning as in Figure 2. 



Figure 5: Pseudorapidity distributions for transverse energy of secondary produced parti- 
cles: all secondaries (Fig. 5a) , all neutrals( Fig-5b) , all charged (Fig. 5c) and gluons (Fig.5d) 
in p + Air Nucleus interactions . The dotted(for 1 TeV laboratory energy) , dashed(for 100 
TeV laboratory energy), and solid(for 1000 TeV laboratory energy) histograms are theoret- 
ical values given by HIJING model. 



Figure 6: Pseudorapidity distributions for p + Air — > tt + + X (Fig. 6a), p + Air — > it + X 
(Fig.6b), p + Air — > K + + X (Fig.6c), p + Air — > K~ + X (Fig.6d) , at 1 TeV,100 TeV and 
1000 TeV. The dotted , dashed and solid histograms have the same meaning as in Figure 5. 



Figure 7: Pseudorapidity distributions for p + Air — >• tt° + X (Fig. 7a), p + Air — >• 7 + X 
(Fig.7b), p + Air -> A + X (Fig.7c), p + Air —> D ± + X (Fig.7d),at 1 TeV,100 TeV and 
1000 TeV. The dotted , dashed and solid histograms have the same meaning as in Figure 5. 



Figure 8: Transverse momenta distributions for p + Air — > all charged + X (Fig.8a), p + 
Air -> p + X (Fig.8b), p + Air -> 7 + X (Fig.8c), p + Air -»• p + X (Fig.8d), at 1 TeV,100 
TeV and 1000 TeV. The dotted , dashed and solid histograms have the same meaning as in 
Figure 5. 

Figure 9: Test of Feynman scaling in the production of p + Air — > 7r ± + X collisions (Fig. 9a) 
and p + Air —>■ p + p + X collisions (Fig.9b), between 1 TeV - 1000 TeV. The Feynman -xp 
distributions were calculated with HIJING model. The dotted , dashed and solid histograms 
have the same meaning as in Figure 5. 

Figure 10: Charged multiplicities distributions in p + Air interactions at 1000 
TeV.Contibutions from all events are depicted in Fig. 10a. In Fig. 10b the histogram is from 
HIJING model calculations with contributions from events with Nj et = 0, Fig. 10c the his- 
togram is from calculations with contributions from events with Nj et = 1 and in Fig.lOd the 
histogram is from calculations with contributions from events with Nj et > 1. 
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Figure 11: Pseudorapidity distributions for all secondaries produced in Fe+Air interac- 
tions(Fig.lla),S+Air interactions(Fig.llb), Ne+Air interactions(Fig.llc)and He+Air inter- 
actions (Fig.lld) at 17.86 TeV/Nucleon for events generated in different impact parameters 
intervals (b min , b max ). See the text for explanations. 

Figure 12: Pseudorapidity distributions for transverse energy of secondary produced particles 
:all secondaries(Fig.l2a); all neutrals(Fig.l2b);all charged(Fig.l2c);gluons(Fig.l2d) in A + 
Air interactions at 17.86 TeV/Nucleon. The theoretical values were calculated with HIJING 
model and are depicted by dotted (He+Air) ; dot-dashed (Ne+Air) ; dashed (S+Air) and 
solid (Fe+Air) histograms. 

Figure 13: Pseudorapidity distributions for transverse momenta of secondary particles for 
A + Air -> ir ± + X (Fig.l3a) and A + Air -> p + X colisions (Fig.l3b) at 17.86 TeV/Nucleon. 
The histograms have the same meaning as in Figure 12. 

Figure 14: Pseudorapidity distributions for A + Air — > ir + + X (Fig. 14a), A + Air — > n~ + X 
(Fig.l4b), A+Air -> K + +X (Fig.l4c), A+Air -> K~+X (Fig.l4d), at 17.86 TeV/Nucleon. 
The histograms have the same meaning as in Figure 12. 

Figure 15: Transverse momenta distributions for A + Air — > all charged + X (Fig. 15a), 
A + Air -> p + X (Fig.l5b), A + Air -> tt+ + X (Fig.l5c), A + Air -> tt" + X (Fig.l5d),at 
17.86 TeV/Nucleon. The histograms have the same meaning as in Figure 12. 

Figure 16: The Feynman xf distributions in the production of A + Air — > 7r ± + X collisions 
(Fig. 16a) and A+Air — > p+p+X colisions(Fig.l6b), at 17.86 TeV/Nucleon. The histograms 
have the same meaning as in Figure 12 . 

Figure 17: Charged particles multiplicities distributions in Fe+Air interactions (Fig. 17a) and 
S+Air interactions (Fig. 17b) at 17.86 TeV/Nucleon. Contributions from events with number 
of minijets Njet > 1 are represented in Fig.l7c for Fe+Air interactions and in Fig.l7d for 
S+Air interactions. 

Figure 18: Pseudorapidity distributions for Fe + Air — > n + + X (Fig. 18a), Fe + Air — > 
7T- + X (Fig.l8b), Fe + Air -> K+ + X (Fig. 18c), Fe + Air -> K~ + X (Fig.l8d) at 
17.86 TeV/Nucleon and impact parameter interval (0-5 fm). Solid histograms are HIJING 
model predictions. Also showing the influence of nuclear shadowing and jet quenching effects 
(dashed histograms). 
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Figure 19: Rapidity distributions for gluons (Fig. 19a) and quark- ant iquark pairs (qq) 
(Fig. 19b) , transverse momenta distributions for valence and sea partons (Fig. 19c) and va- 
lence parton (Fig.l9d) for Fe+Air interactions at 17.86 TeV/Nucleon and impact parameter 
interval (0-5 fm).The solid and dashed histograms have the same meaning as in Figure 18. 



Table 1: Particle multiplicities for pp interaction at 200 GeV are compared with data from 
Gazdzicki and Hansen |32|] and with the results given by VENUS model (as computed here) 
and DPM model DPMJET II 0. 



PP 


Exp. data 


HIJING 


HIJING U) 


VENUS 


DPMJET II 


< 7T~ > 


2.62 ±0.06 


2.61 


2.65 


2.60 


2.56 


< 7T+ > 


3.22 ±0.12 


3.18 


3.23 


3.10 


3.17 


<7T° > 


3.34 ±0.24 


3.27 


3.27 


3.28 


3.38 


< h~ > 


2.86 ±0.05 


2.99 


3.03 


3.05 


2.82 


<K + > 


0.28 ±0.06 


0.32 


0.32 


0.27 


0.28 


< K~ > 


0.18 ±0.05 


0.24 


0.25 


0.19 


0.19 


< A + S u > 


0.096 ±0.015 


0.16 


0.165 


0.18 




< A + S° > 


0.013 ±0.01 


0.03 


0.037 


0.033 




<K°> 


0.17 ±0.01 


0.26 


0.27 


0.27 


0.22 


<p> 


1.34 ±0.15 


1.43 


1.45 


1.35 


1.34 


< p > 


0.05 ±0.02 


0.11 


0.12 


0.06 


0.07 



Table 2: Particle composition of p + p interactions at 540 GeV in cms. 



Particle type 


< n > 


Exp. data 


HIJING^ 


All charged 


29.4 ±0.3 





28.2 


K°±K° 


2.24 ±0.16 


in 


1.98 


K+ ±K 


2.24 ±0.16 


H 


2.06 


P±P 


1.45 ±0.15 


@ 


1.55 


A ± A 


0.53 ±0.11 




m 


0.50 


£+±£-±E+±E- 


0.27 ±0.06 




m 


0.23 


h- 1 


0.04 ±0.01 




|3j 


0.037 


7 


33 ±3 


in 


29.02 


7T + ± 7T~ 


23.9 ±0.4 





23.29 


K° 

s 


1.1 ±0.1 


H 


0.99 


7T° 


11.0 ±0.4 


@ 


13.36 
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Table 3: Percentage of occurence of various particles among the secondaries of a Nucleus- 
Air collision as calculeted by the HIJING and VENUS models. The number of protons and 
neutrons have been reduced by respective numbers in the primary system(values labeled by 
starf*) 



201). The average multiplicity and numbers of participants are also given 



Particle 


Projectile 


56 Fe 


aa S 


™Ne 


4 He 


P 


P 


type 


Energy (TeV/N) 


17.86 


17.86 


17.86 


17.86 


17.86 


1000 


< 7T~ + 7T + > 


HIJING 


45.98 


45.93 


45.96 


45.86 


45.76 


46.57 




VENUS 


51.02 


51.29 


51.48 


52.34 


53.05 


52.15 


<7T° > 


HIJING 


26.13 


26.02 


26.10 


25.93 


26.04 


26.38 




VENUS 


28.30 


28.52 


28.49 


28.85 


28.93 


28.43 


< K+ + K :> 


HIJING 


5.20 


5.20 


5.20 


5.20 


5.16 


5.58 


K mesons 


VENUS 


12.35 


12.02 


11.84 


11.11 


10.51 


10.86 


<K°> 


HIJING 


2.54 


2.57 


2.53 


2.52 


2.57 


2.75 




VENUS 














< P > 


HIJING 


4.16 


4.20 


4.20 


4.37 


4.78 


3.25 


< p > * 


VENUS 


0.70 


0.66 


0.66 


0.66 


-0.11 


0.66 


< n > 


HIJING 


4.12 


4.20 


4.23 


4.40 


4.16 


2.86 


< n > * 


VENUS 


0.54 


0.61 


0.64 


0.66 


1.49 


1.43 


_ A _ 

< A > 


HIJING 


0.70 


0.70 


0.71 


0.72 


0.72 


0.56 


<. iv ~r Z-i > 




1.00 


1.041 


1 fi7 
1.0 ( 


1 fin 

1.0U 


1.00 


l.Ol 


Other baryons 


TTT TTAT/~1 

H1J1JNG 
















VFNTT^l 

V H/1M KJ o 


n zLfi 
u.^to 




U.^tZ 


n ^zl 


n 9fi 
u.zo 


n 98 
u.zo 


< 7 > 


TTT TTM/^ 


a on 

4.30 


a on 

4. JO 


A O O 

4.33 


A O O 


A O *T 

4.27 


A AO 

4.48 




V H/1M KJ o 


1 fin 

l.UU 




1 A.'X 
1 .'-to 


1 90 
l.ZU 


1 1 9 


IIS 
i.ia 


111 J 

all charged 


TTT TTAT/~1 

HIJING 


f T O n 

57.30 


r - T A T 

57.47 


57.54 


r' T f A 

57.54 


f *1 O A 

57.84 


57.79 


all neutrals 


■ ■I 1 1 rv f 1 
111,1 liN \j 


4:2.66 


4z.e> ( 


4z.ol 


4z.oU 


a o /in 
4z.4U 


/in i n 

4z. lo 


< P > 


HIJING 


1.42 


1.38 


1.43 


1.40 


1.41 


1.56 


< n > 


HIJING 


1.40 


1.38 


1.43 


1.40 


1.41 


1.60 


< gluons > 


HIJING 


9.70 


9.20 


8.96 


7.73 


7.20 


10.71 


< q + q > 


HIJING 


0.40 


0.43 


0.41 


0.36 


0.33 


0.71 


Mean 


HIJING 


274.0 


217.5 


203.0 


88.7 


38.9 


51.7 


multiplicity 


VENUS 


354.7 


270.5 


209.3 


92.6 


49.4 


106.5 


Mean projectile 


HIJING 


9.6 


7.2 


6.5 


2.2 


1.0 


1.0 


participants 


VENUS 


12.1 


8.1 


5.7 


2.1 


1.0 


1.0 


Mean target 


HIJING 


5.9 


5.6 


5.6 


3.4 


2.06 


1.42 


participants 


VENUS 


6.0 


5.3 


4.7 


3.2 


2.0 


2.06 
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Table 4: Mean transverse energy for the secondaries of a Nucleus-Air collisions at 17.86 
TeV/Nucleon as calculeted by the HIJING model and by considering superposition of N co u 
nucleon-nucleus collisions and N pro j nucleon-nucleus collisions,where N co u is the number of 
binary collisions and N pro j is the number of participant projectile nucleons. E^ A is mean 
transverse energy in p+Air interaction at the same energy and in the same impact parameter 
interval (see the text for explanation). 





Projectile 


56 Fe 


32 S 


20 Ne 


4 He 


P 


Mean number 


< N coll > 


27.75 


19.21 


13.61 


2.62 




Mean number 




18.43 


11.44 


8.24 


1.73 




Mean Transverse 


all secondaries 


212.5 


140.7 


104.2 


25.16 


7.54 


Energy 


Ncou * E£ A 


209.2 


144.8 


102.6 


19.75 




(GeV) 


Nproj * E T 


138.96 


86.25 


62.10 


13.04 




Mean Transverse 


all charged 


122.7 


81.27 


60.10 


14.52 


4.39 


Energy 


N C oii * E£ A 


121.8 


84.33 


59.74 


11.50 




(GeV) 


Nproj * E^ 


80.90 


50.22 


36.17 


7.6 




Mean Transverse 


all neutrals 


89.87 


59.43 


44.10 


10.64 


3.15 


Energy 


Ncoii * EP A 


87.41 


60.51 


42.87 


8.25 




(GeV) 


Nproj * Erjr 


58.05 


36.03 


25.95 


5.44 
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P + P(300 GeV) RAPIDITY AND TRANSVERSE MOMENTUM H II J II NG 
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